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Abstract. This paper deals with a half-space linearized problem for the distribution function of 
the excitations in a Bose gas close to equilibrium . Existence and uniqueness of the solution, as well 
as its asymptotic properties are proven for a given energy flow. The problem differs from the ones 
for the classical Boltzmann and related equations, where the hydrodynamic mass flow along the 
half-line is constant. Here it is no more constant. Instead we use the energy flow which is constant, 
but no more hydrodynamic. 

S3 ■ 1 Introduction. 
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in 



We consider a linearized half-line problem related to the kinetic equation for a gas of excitations 
interacting with a Bose condensate. Below the temperature T c where Bose-Einstein condensation 
sets in, the system consists of a condensate and excitations. The condensate is modelled by a 



Gross-Pitaevskii equation. The excitations are described by a kinetic equation with a source term 
taking into account their interactions with the condensate, i.e. 



O ■ at 
m 



dF 

+ p . VxF = Cl2 (F,n c ). (1.1) 



With F the distribution function of the excitations, and n c the density of the condensate, the 
collision operator in this model is 

C 12 (F,n c )(p) = n c J x$ c 5 3 ((1 + F 1 )F 2 F 3 - F 1 (l + F 2 )(l + F 3 ))dp x dp 2 dp^ (1.2) 

where 

5 C = 8 (pi =P2 +P3,Pi = p\ +P3 + n c), h = 8(p = pi) - S(p = p 2 ) - 8(p = p 3 ), 

and x is the characteristic function of the set of (p,pi,P2,P3) such that for some A 3> n c , 

\p\ < A, |px| < A, < A, |p3| < A. 

For a discussion of this two-component model, we refer to pQ and references therein. The kernel of 
C\ 2 consists of the Planckian distribution functions 

P a b(p) = ? 2 ^ \- , P G M 3 , for a > 0, /3 € K 3 . 
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This paper considers a linearized half-space problem connected to Assuming F near the 

Planckian equilibrium 

Pc(p) ' 



e p 2 +n c _ I ' 



write F = P c (l + /). For simplicity we denote n c (resp. P c ) by n (resp. P) below. This gives the 
integrand of the collision operator 

(1 + F l )F 2 F 3 -F 1 (l + F 2 )(l + F 3 ) = - (1 + P 2 + P 3 )P 1 f 1 + (P 3 - Pi)P 2 / 2 + (P> - Pi)P 3 h 

+ P2P3/2/3 — -P1-F2/1/2 — P\P$hh- 

Here we have used that (1 + P) = M~ 1 P, where 

M(p) = e" p2 " n , pes 3 , 

and that M(px) = -^(p 2 )-^(p3) when p\ = p^+p^ + n holds. It follows that the linearized term in 
the previous integrand gives the linearized operator 



Tt f 

L(f) = p J XS C S 3 I - (1 + P2 + P 3 )Pi/i + (Ps ~ P1W2 + (P2 - Pi)P 3 f 3 



dpidp 2 dp 3 . 



We consider a half-line in the x-direction, and functions which are cylindrically symmetric in the 
p- variable, i.e. functions of p x and with only radial dependence p r in the orthogonal direction. This 
changes the momentum conservation Dirac measure in C\ 2 and L to 5(p\ x — p 2x —p 3x )- 
The Milne problem is 

Pxd x f = Lf, x> 0, p x G R, Pr e M + , (1.3) 



/(0,p) = /o(p), Px>0, (1.4) 

where fo is given. We will prove in Section 2 that KerL is spanned by (p 2 + n)(l + P) and p x (l + P). 
For any measurable function f(x,p) such that for almost all x E M + , 



H/(^))6ip (1+|p|) a( 



let 
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/(a:,p) = a(x)(p 2 + n)(l + P) + b(x)p x (l + P) + w(x,p), (1.5) 
be its orthogonal decomposition on the kernel of L and its orthogonal complement in L 2 P , i.e. 

T+P 

J p x w(x,p)Pdp = J (p 2 + n)w(x,p)Pdp = 0, x £ R + . (1-6) 
The main result of this paper is the following. 
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Theorem 1.1 For any £ G R and f G L 2 (1+ | p | )3 _^ (R + x R+), there is a unique solution f G D 
to the Milne problem, 

Pxd x f = Lf, x> 0, Px G R, Pr G R+, (1.7) 

/(0,p) = /o(p), Px>0, (1.8) 

p x {p 2 + n)f(x,p)P(p)dp = £, x€R + . (1.9) 
Moreover, for the decomposition of the solution, there are (ooo^oo) £ I^ 2 suc/i £/ia£ 

6oo = -, (1.10) 

7 

and a constant c > 0, sac/i i/iat /or any n G]0, C4[, 

(l + |p|) 3 n; 2 (x,p)^pdp + |a(x)-a 00 | 2 + |6(2;)-6 00 | 2 <ce- 2r ' x , x G R + . (1.11) 

Here, Z? is the function space defined by 

D = {f;fe L 00 (R+;Lj ) (1+|p|)3 _p_(R x R+)) 5 p x d x f G lL(R + ! ^ r(1+|p|) - 3 _p_ (K X R + ))}, 

7 = j P Up 2 + n)P(l + P)dp, c 3 = ^j piP(l + P)dpj pl(p 2 + n) 2 P(l + P)d P y, 
c 4 = mzn{y,— }. (1.12) 
Remarks. 

ITus result should be compared to the analogous result concerning the Milne problem for the lin- 
earized Boltzmann operator around the absolute Maxwellian in In ^ the mass flow is constant 
and well-posedness for the Milne problem is proven for a given mass flow. In the present paper on 
the other hand, the mass flow is no more constant, since mass is mot a hydrodynamic mode. But 
the energy flow is constant, and well-posedness here is proven for a fixed energy flow. That this 
energy flow is proportional to the asymptotic limit of the mass flow is a new low temperature result. 
A separate complication in the present case is that, whereas the given mass flow in U\[ is a hydro- 
dynamic component of the solution, this is no more the case for our given energy flow. Another 
differing aspect compared to classical kinetic theory, is that the collision frequency is asymptotically 
equivalent to \p\ 3 , when p — > oo. 

The interest in half space problems such as (jl.7p - (|1.8p is partly due to their role in the asymptotic 
behaviour of the solution of boundary- value problems of kinetic equations for small Knudsen num- 
bers. This subject has received much attention for the Boltzmann equation |12| |14| |15j |16] |17j 
[2] and related equations [6] [13]. Starting from the stationary Boltzmann equation in a half-space 
with given indatum and a Maxwellian limit at infinity, the unknown is assumed to stay close to 
this Maxwellian, giving rise to the linearized stationary Boltzmann equation in a half space. The 
general treatment of the linearized problem for hard forces and hard spheres under null bulk ve- 
locity, is given in |llj and references therein. The case of a gas of hard spheres (resp. of hard or 
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soft forces) and a null bulk velocity at infinity is independently treated in [2] (resp. in [IQ]). The 
case of a gas of hard spheres and a nonzero bulk velocity at infinity is considered in [8j, positively 
answering a former conjecture [7]. The Milne problem for the Boltzmann equation with a force 
term is analyzed in |9j. Half-space problems in a discrete velocity frame are studied in [3]. For a 
review of mathematical results on the half-space problem for the linear and nonlinear Boltzmann 
equations, we refer to [5]. 

The plan of the paper is the following. In Section 2, the linearized collision operator L is stud- 
ied, including a spectral inequality. In Section 3, Theorem II .11 is proven. 

2 The linearized collision operator. 

Lemma 2.1 

L is a self-adjoint operator in L 2 P . Within the space of rotationally invariant distribution func- 

T+P 

Hons, its kernel is the subspace spanned by (p 2 + n)(l + P) and p x (l + P)- 
Proof. It follows from the equalities 

P 2 (l + P 2 )(P 3 ~ Pi) = P 3 (l + P3)(P2 ~ Pi) =Pl(l + P 2 )(l + Ps), 

p(l,p>p\-PP- P l( 1 + - P 2)(l + P3) 2 _ _2 , _2 , _ 

Pl{L + f 2 + f 3 ) — f 2 f 3 — -i , 5 1 Pi — £>2 + P3 + n > 

1 + f\ 



that for any functions / and g in L 



2 

p i 
i+p 



J Y^(P)f(v)Lg(p)dp = -n J x<5 c Pi(l + P 2 )(l + P 3 



fi / 2 h \ f 91 92 93 A , , , 

-)( -, , „ - -i , „ - ; , „ )dp x dp 2 dp 3 . 



V l + Pl 1 + P 2 1 + P 3 /V 1 + P! 1 + P 2 1 + P 3 

This proves the self-adjointness of L in L 2 P . Moreover, Lf = for / G L 2 P implies that 

T+p T+P 

i T~p = 1 + I ' Pi* =P2x +P3x, P?=P2+P3+ n - 

1 + -Tl 1 + -T 2 1+^3 

It is a consequence of this Cauchy equation that the orthogonal functions 

(p 2 + n){l + P) and p x (l + P) (2.1) 
span the kernel of L. 



The operator L splits into K — v, where 

Kf(p) := p£y( J xKPx = P2x+P3x,P 2 =P2+P3 + n )( p 3 - P)P2f2dp 2 dp 3 
+ J X$(pix = Px + P3x,Pi = P 2 + p\ + n )(! + P + P 3 )Pifidpidp 3 

+ y X*(Pl* = + Ps^P 2 = P 2 + p\ + n)(Pi - P)P 3 f 3 d Pl dp 3 ^j (2.2) 



and 



v(p) := n J x5(p x = p 2x + P3x,p 2 = pI+pI + n)(l + P 2 + P 3 )dp 2 dp 3 

+2n J x${Pix = Px +P3x,p( =P 2 +pI + n)(P 3 - P{)dpidp 3 . (2.3) 

Notice that K = K\ as well as v = v\ depend on A through % 
Lemma 2.2 

K\ is a compact operator in L 2 P . Moreover, there are positive constants vq and v\ , such that 

T+P 

uniformy with respect to A, 

^o(l + H) 3 <v\(p) <^i(l + H) 3 , P =(fc,p r )eKxR+, for |p| < A. (2.4) 

Proof of Lemma 2.2. 

K = K\ + K 2 + K 3 , where 



n) — - — dp 3 , 



#i/(p) := 2n J X 1 (p,p 2 )P 2 f 2 dp 2 , 
Xi(p,p 2 ) := J x5{p x =p 2x +p 3x ,p 2 =pl+pl + 
K 2 f(p) := 2n j X 2 (p,p 1 )P 1 f 1 dp 1 , 

Xd~{pix =Px+P3x,pj =p 2 +pl+n) -dp 3 

K 3 f(p) := 2n J X 3 (p,p 3 )P 3 f 3 dp 3 , 

/p p 
X$(p lx =p x +p 3x ,pj =p 2 +pl + n)-^ — dpx. 



Then 



Xi{p,p 2 ) = p^(P(\Jp 2 -P 2 ~n) - P{p))Xp2-pl- n >0 j XXpl> P l+pl x +J{Px = P2x + P3x)dp 3x 

7T 1 _ 

= Pjp)\p 2 -P% _ 1 " F W» XX ^P2x-b*) 2 +Plr<lvl+vl-ni 

7T 1 
X 2(P,P1) = p^j(l + ^(P) + ep 2_ p2 _ 1 )XX p j-pg-( pig - BB )g-n>0. 

7T 1 

We see from the above that K\ (resp. K 2 , K 3 ) is bounded in I? P if and only if A is finite. When 

T+P 

the domain of p is bounded, K\ is compact as an Hilbert-Schmidt operator. K 2 and K 3 behave 
similarly, and so if is a compact operator in L 2 



p ■ 

i+p 



Since 

P 3 - P 1 = Pi(l + P 3 )(e p2+n - 1) > Pm, p? = p 2 + pi + n, 
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the function v is such that 



v(p) > 2n 2 j X\ pi \<\,\ P3 \<\ p i${Pix =Px+P3x,pI =P 2 +pl + n)dp 1 dp 3 
= 2vrn 2 / X\ Pl \<\Pidpi 

J p\~P 2 -(Plx-Px) 2 -n>0 

= 2vrn 2 / X\ P1 \<\Pidpi 

J p\ r +2p x p lx >2pl+pl+n 

> 2vrn 2 / X\ P1 \<\Pidpi, p 2 < 100n. 

Jv<rVi<r>0. p?„>201n 



'PxPix>0, pf r >201n 

Moreover, for A 2 > p 2 > lOOn, the first term of is equal to 

J X\p 2 \<\\pz\<\KPx = P2x+P3x,p 2 =pI+pI + n)(l + P2 + Ps)dp2dp 3 

r , r\/^ 2 -( P 2x- P x) 2 

> 4VT 2 / ( / (5(p§ r =p 2 r -p\- P\ x + 2 P*P2z - n)p 3r dp 3r )p2rdp2rdp2x 

J\p2x—Px\<\P2r>0,\P^\<^ J 

= 27T 2 / P2rdp2rdp2x 
J\P2x~Px\<\ P2r>0, |P2|<A, 2p| a .+p| r -2p a; p2 :c <P?-n 

> 27T 2 / P2rdp2rdp2x 

J\P2x-Px\<\ P2r>0, PxP2x>0, 2p| a .+p| r .+n<p2+2p, c p2 ;c , n<p^<4n 

• > 27T 2 / P2rdp2rdp2x 

Jp2r>0, PxP2x>0, 9n<pl+2p x p2x, n<p? 2 <An 

>2TT 2 P2rdp2rdp2x 
Jp2r>0, PxP2x>0, 9n<p2+2p x p2x , n<p\<kn, p 2 x >n 

>2TT 2 P2rdp2rdp2x- 
Jp2r>0, pxP2x>0, n<p?,<4n, p\ x >n 

And so, v is bounded from below by a positive constant. 
The first term of u(p) belongs to the interval with end points 

2vr 2 nx| p |<A / _ P2rdp2rdp2x 
and 



P2r>0,pl r +2{p2x-\px) 2 <\p%+p 2 --n,\p2x~Px\<WP2\<\ 



27T 2 n(l + — )X|p|<A / P2rdp2rdp 2 x- 

e 1 Jp2r>0,p 2 2r +2{p2x-\px) 2 <\p 2 +p 2 -n,\ P 2x-px\<X\p2\<\ 
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Then for \p\ < A , with the change of variables (x,y) := {j>2x->P2r)-> 



2 / P2rdp2rdp2x 
Jv^>0,v 2 ^+2(v^-\p x ) 2 <\% 



lp2r>0,p^ r +2(p 2x -^p x ) 2 <^p 2 +p 2 -n,\p 2x -p x \<X,\p 2 \<\ 



= / dxdy 
Jy^ 



!y>0,y+2(x-±p x ) 2 <±p 2 +p 2 -n,\x~p x \<\,x 2 +y<\ 2 



, 1 '2 , _2 



r2P x +Pr-n r 

Jo J( x -4)p x ) 2 <h 



dxdy 



-L 



yjp 2 x + 2p2 -In- 2ydy 



= \{pl + 2pl-2n) 3 2 

~ H 3 . 

It has here been used that for any positive A and any \p\ < A, 

y + 2(x - ^p x ) 2 < -pi +p 2 r =► |x -p x \ < A and x 2 + y<A 2 . 
The second term of z/(p) is bounded. Indeed, for \p\ < A, 

< J X\ Pl \<\,\p 3 \<\$(Pix = Px+P3x,Pi =P 2 +pl + n)(P 3 - P{)dpidp2, 

< [ X| P i|<a( ~2 i 7 - A) ( / = Pir - - + 2 PxPis - ri)p 3r dp 3r ) p lr dp lr dp lx 

J Plr >o e p i p - 1 v Jo 7 

- 4 J J eX 2+y- p 2 _ ^Xy>p2+pl+n-2xp x dydx 

= \Y, ekl>2 I ^ l e ~ ky Xy>pl+pl+n-2xp x dydx 
k>l J J ° 

i J 



k>i 



k>l 



Denote by (-,-) the scalar product in L 2 _^, and by P the orthogonal projection on the kernel 

T+p 

of L. 

Lemma 2.3 

L satisfies the spectral inequality, 

~(Lf, f) > H)((l + H) 3 (/ - P)f, (I - P)f), f G 4 +W) 3_p_- (2-5) 

Proof. For the compact, selfadjoint operator K, the spectral inequality holds as in the corresponding 
Boltzmann case. Namely, there is no eigenvalue a > 1 for ^. Else there is / / such that 
Lf = (a- l)vf and so (Lf, f) > 0. But 



(LfJ) = -nj 5 C (-J^ - -J^- - -^-?d Pl dp2dp, < 0. 
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Except for the kernel of L discussed above, the eigenvalues of are bounded from above by «o < 1- 
Spanning L 2 3 with the eigenfunctions of we obtain the spectral inequality 



(Lf, f) < (a - l)W - P)f, (I - P)f), f € L* 



(i+H) 



3 P • 
l+P 



Finally, ([23]) follows from (J2 

3 The Milne problem 

In this section, we prove Theorem ll.il 
Let 

/ = / - -Px(i + p), Hp) = Hp) - -Px(i + p). 

7 7 
Solving the Milne problem (jl.7jl - (|1.8p - (|1.9p for the unknown / is equivalent to solving 

PxdJ = Lf, x > 0, V x eR, Pr €R + , (3.1) 
f(0,p) = Hp), Px>0, (3.2) 

J p x (p 2 + n)f(x,p)P(p)dp = 0, x£R + , (3.3) 
for the unknown /. 



We study the behaviour of the solution / to the Milne problem (|3.1]) - (|3.2p - (|3.3p when x — > +oo. 
Set 

f{x,p) = (a(x)(p 2 + n) + b(x)p x )(l + P) + w(x,p), 

with 

J p x wPdp = J (p 2 + n)wPdp = 0, 
an orthogonal decomposition of /. Denote by 

W(x) = 1 J Px p(x,p) T ^dp (3.4) 

the linearized entropy flux of /. It holds that 

W(0) < l -f Px f 2 (p)-J—dp. (3.5) 

By (3.3) 

W{x) = \ J PxP(x,p) 1 ^ p dp- ^ J plf(x,p)P(p)dp J p x (p 2 + n)f(x,p)P(p)dp 
= 2 / P xVj2 ( x >p) i _[_ + a y Px(p 2 + n )wPdp + b J p 2 wPdp + ab^ 
— ^ ( y p 2 x wPdp + aj^j ^ y pa (p 2 + n)wPdp + 67^ , 



i.e. 



W(x) = — J p x w 2 (x,p) - — p dp j p 2 x wPdp j p x (p 2 + n)wPdp. (3.6) 

This differs from ([BCN]), where the linearized entropy flux of the solution is equal to the linearized 
entropy flux of its non-hydrodynamic component. The expression (3.6) of W in terms of w is impor- 
tant for the proof. 

Multiplying (|3.ip by fjrj>, integrating on (0, X) x R x R + (resp. R x R+) and using (|2.5p . gives 
W(X) + v J(l + \p\fw 2 (x,p) T ^dpdx< W(0), X>0, (3.7) 

and 

W\x) + u Q J(l + \p\f w \ X) p)-JL_dp < 0. (3.8) 

Since / € D, it holds that W G L°°(R+). Then by ([321) and W G ^ 1 ( R+ )- B y *M ? it is a 

non increasing function. Hence it tends to zero when x tends to +oo and is a nonnegative function. 
Let r] S]0,c 4 [. Multiply (33} by e 2r?:r , so that 

{W{x)e 2r > x )' - 2rjW(x)e 2 ^ x + i^e 2 ^ j (I + |p|) 3 w 2 (x,p) y^pdp ^ °- 
By the Cauchy- Schwartz inequality, 

I j p 2 x w(x,p)Pdp J p x (p 2 + n)w(x,p)Pdp\ < ^ J w 2 (x,p) ^ p dp. 
Hence, 

(W(x)e 2 ^)' + e 2r > x J (Ml + \p\) 3 -v(Pz + c 3 ))w 2 (x,p) T ^dp<0, x>0. (3.9) 

By the definition (|1.12p of C4, the nonnegativity of W and (I3.5p . it holds that 

\ e 2r > x / (1 + \p\fw 2 (x,p) -dpdx<c, (3.10) 

Jo J 1 + P 

for some constant c. Moreover, by (|3.9|) and (|3.8p . 

< W(a:) < W(0)e- 2r?a; < ce" 2 ^, x > 0. (3.11) 

(|3.10p implies that f(x, •) converges to a hydrodynamic state when x — > +00. In order to prove 
the exponential pointwise decay of J(l + \p\) 3 w 2 (x,p) j^rpdp in (|l.lip . let < Y < X be given and 
introduce a smooth cutoff function <£(x) such that 

= 0, x e [0, 2"[U]X + l,+oo[, $(a;) = 1, a? € [Y,X]. 
Denote by </?(x) = e r?x $(x). Then, 

Pa,^/) = L(d x (<pf)) + ip'Lw+p xV "f. (3.12) 
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Multiply (13TT2]) by d x (ipf)j£p, integrate over M Pa . x M+ r and use (|2.5p . Hence, 

55£ / P^ d ^f)) 2 TTp dp + U0 J ^ + \p\) 3 (9 x (^)) 2 T ^pdp 



< J d x (<pf)^p'Lw +p x <p"f) - ^ p dp 



ip' j d x ((pw)Lwj^-pdp + <p<ff' j Pxfd x f j ^p d P + v'v" j PxP ^ p dp, 



i.e. 



^ j p^M)) 2 ^dp + v J {i + \p\)%d x ^)f^dp 

< <p' J d x {( P w)Lw^- p dp+ WW)' + w - W (3) )^. 
Integrate the last inequality on [0, +oo[, so that 







p 



//„ / l{l + \p\f{d x {ipw)) 2 ^-^dpdx 



p 



+00 



< / ip'(x) d x (tpw)Lw -dpdx+ / {ip'ip" - ipip (3) )W(x)dx 

Jo J 1 + P Jo 

p+oc r p 1 r+oa r -i p 

r+00 

+ / OV - yyW)W{x)dx 



Jo 

p+oo r p p+oa i- p 

<n <j to f J y a + MmM) 2 — dpd* + — y y (1 + ^v*—^ 
/•+00 

+ / - w (3) )W(s)dx, a>0. 

JO 

Choose a < ,, . Use (|3.10|) . the Cauchy- Schwartz inequality and the exponential decay of W 

jyP Hoc 

expressed in (|3.1ip in the W-term. It then holds 



-OG 







(1 + \p\f(d x {i P w)) 2 -^ p dpdx < c, 



for some positive constant c. Finally, 

e 2vX J (l + \p\fw 2 (X,p) T ^dp = 2j*J(l + \p\f(d x &w)) 2 T ?- p -dpdx <c, X>0. (3.13) 

The exponential decay of (a, b) to some limit (cioo, froo) when x tends to +00 can be proved as follows. 
The solution f(x,p) = (a(x)(p 2 + n) + b(x)p x )(l + P) + w(x,p) is solution to (|3.ip iif 

(a'p^ 2 + n) + Vp 2 x ){l + P) +p x d x w = Lw. 

Multiply the former equation by p x P (resp. (p 2 + n)P) and integrate with respect to p, so that 

(a + — p x w(-,p)Pdp)' = (b + - p x {p 2 + n)w(-,p)Pdp)' = 0. 
7i 7i 
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Denote by 

Coo := o(0) + ~ y p 2 x w(0,p)Pdp, ■= 6(0) + ~ y PxO 2 + n)u;(0,p)Pdp. 

By the Cauchy-Schwartz inequality and fj3 . 13f) . 

i 

|a(x) - aoo| = ~| / p 2 w(x,p)Pdp\ <c(f(l + \p\) 3 w 2 (x,p)j^pdp) 2 < ce~ r,x , 



|6(x)-6oo| = ~| y ^(p 2 + n) W (x,p)P(ip| <c( J(l + \p\) 3 w 2 (x,p)j^d P y <ce^ x . (3.14) 
By HUB 

lim / p x f 2 (x,p) P dp = 0. (3.15) 
By iH) 

lim [(l + \p\) 3 w 2 {x,p)—^—dp = 0. (3.16) 

z-s>+oo y i + p 

Using the decomposition / = (a(p 2 + n) + bp x )(l + P) + w of / into its hydrodynamic and non 
hydrodynamic components and 

doc = lim a(x), boo = hm 6(x), 

>+oo X — ¥+oo 

it follows from (l3TT5]) - (l3TT6l ) that 

lim [ p x ((a(x)(p 2 + n) +b(x)p x )(l + P)) - dp = 0, i.e. 000600 = 0. 
>+oo J \ / 1 -\- Jr 



Below this will be improved to (1-10) 600 = 0. But first we prove the existence of a solution / E D to 
the Milne problem ()3.ip - (|3.2p - (|3.3p . It will be obtained as the limit when / — > +00 of the sequence 
(fi)i€N* of solutions to the stationary linearized equation on the slab [0, /] with specular reflection 
at x = I, i.e. 

PxdJi = Lf u x e [0, 1], Px e R, Pr e R+, (3.17) 
fi(0,p) = Up), p x >0, (3.18) 

fl(l,Px,Pr) = fl(l,-Px,Pr), Px < 0. (3.19) 

We switch from given indata and no inhomogeneous term to indata zero and an inhomogeneous 
term. Let e > be given. Let the subspace D(A) of L 2 „ P ((0, /) x R x R+) be defined by 

Pr(l + \p\) jZfp 



D(A) = {ge L 2 pril+lp{)3 ^{{0,l) x R x M+); Px d x g E ^ r(1+w)3 ^ ((0, 1) x R x R+), 



^ + ); rA<? E L 2 

1+P Pr\^\i>\j 1+p 

g(0,p) = 0,p x > 0, g(l,p x ,Pr) = 9(1, ~Px,Pr),Px < 0}. 
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The operator A defined on D{A) by 

(Ag)(x,p) = eg(x,p) +p x d x g(x,p) 

is m-accretive since I — ^A is bijective. Indeed, for any /6lJ (i+| p |)3 p ((^' x ^ x ' there 
is a unique g £ D(A) such that 

(I-±A)g = f i.e. ^g-^Pxd x g = f. (3.20) 



Here 5 is explicitly given by 

2e P c ^-y 

s(x,p) = / f(y,p)e dy, p x > 0, 

Px Jo 
rl 



2e / f x+y-21 r l x- y \ 

g(x,P) = —y f(y,p)e e dy + / f(y,p)e t p^dyj, p x <0. 



It belongs to L 2 (1+ | p | )3 _p_(( >0 x K x M+ ) since multiplying (^20]) by 2p(l + IpI) 3 ^, then 
integrating on [0,1] X M 3 implies that 



/ g 2 {x,p){l + \p\f-^—dxdp + ^- f \p x \(l + \p\fg 2 (0,p)-—-d P 

J ' Ze Jp x <0 1 + -r 

= 2 J f( XlP )g( X ,p)(l + \p\f-JL_dp 



< J f 2 {x,p){l + \p\f-^pdxdp + j g 2 (x,p)(l + \p\) 3 -^-pdxd P . 

It then follows from (pT20l) that p x d x g 6 L 2 . P ((0, i)xRx 

— L is an accretive operator. Hence A — L is m-accretive and there is a unique solution 

/ e €Lj p(1+|p|) 3 3 p p ((0,OxRxR+) to 

e/e + = L/ e , x > 0, P:c 6 K, p r 6 R+ (3.21) 

/ e (0,p) = /o(p), Px>0, 

fe(l,Px,Pr) = fe(l,-Px,Pr), Px < 0. 

In order to prove that there is a converging subsequence of (/ e ) when e tends to zero, split f e into 
its hydrodynamic and non-hydrodynamic parts as 

f e (x,p) = (a e (x)(p 2 + n) + 6 e (x)p a; )(l + P) + w e (x,p), 

with 

y p x w £ Pdp = J(p 2 + n)w e Pdp = 0. (3.22) 

Multiply (j!T2T]) by / eI ^p, integrate w.r.t. (x,p) G [0, I] xMx R + and use the spectral inequality (|2.5|) , 
so that (w e ) is uniformly bounded in L 2 ( 1+ [ p Q3 p ^] x R x R + )- Notice that the boundary term 

at 2 vanishes. And so, up to a subsequence, (w e ) weakly converges in 

l2 ,ui ns p ([0,Z]xRxR+)) 
to some function to. Moreover, the same argument as for getting ()3.13j) can be used here, so that 

e 2vx J(l + \p\) 3 w 2 (x,p) T ^dp<c, xe[0,l}. (3.23) 
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Expressing f px>0 p x fe{0,p)j^pdp (resp. f Px>0 Px{p 2 + n)f € (0,p)j^pdp) in terms of o e (0), & e (0) and 
w e {0, •) leads to 



and 



a e (0) / Px(p 2 + n)Pdp + 6 e (0) / p 2 x Pdp 

= I Pxfoip) --. P v d P- [ PxW e (0,p) P dp, 

Jp x >0 1 + -r J Px >o 1 + r 

o e (0) / p x (p 2 + n) 2 Pdp + 6 e (0) / p 2 x {p 2 + n)Pdp 

= [ Px(p 2 + n)f (p) P dp- f p x {p 2 + n)w e {0,p) P dp. 
Jp x >0 1 + -r J Px >0 x + Jr 

By the Cauchy- Schwartz inequality and (|3,23p taken at x = 0, it follows that 
/ p x u» e (0,p) P dp and / p z (p 2 + n)w e (0,p) P dp 

are bounded. Consequently, (a e (0), 6 e (0)) is uniformly bounded with respect to e. Moreover, f e 
solves (pmD iif 

e((a e (p 2 + n) + iy^Xl + P) + w e ) + «Px(p 2 + n) + 6^ 2 )(1 + P) +p a .0 !B ii> e = Lu» e . (3.24) 
Multiplying the previous equation by p^P (resp. (p 2 + n)P) and integrating w.r.t. p implies that 

eb e J p 2 x P(l + P)dp + We + ( / plwePdp)' = 0, 
ea t J (p 2 + n) 2 P(l + P)dp + -yb'e+(j Px (p 2 + n)w e Pdp}' = 0. 
Consequently, denoting by 



a 



' j p|P(l + P)dp and (3 = J J {p 1 + n) 2 P(l + P)dp, 



it holds that 



7a e (x) = — J p 2 .w € {x,p)Pdp + (7a £ (0) + J p x w e {0,p)Pdp)e t x 

+ p x ?—(p 2 + n)w t {y,p)Pdp^e~ ( - x ~ y ^dy, x£ [0,1], 

'jb e ( x ) = ~ J Px(p 2 + n)w e {x,p)Pdp + (7& e (0) + / p x {p 2 + n)w t {0,p)Pdp)e" 
( y — P 2 w e (y,p)Pdp^e ~^ x ~ y ^dy, xG [0,1]. 



Together with the bounds of (a e (0), 6 e (0)) and (|3.23[) . this implies that (a e ) (resp. (b e )) is bounded 
in L 2 . And so, up to a subseque: 

solution /, of (|3T1T1) - (I3TT8^ - ([3TT91) . 



in L 2 . And so, up to a subsequence, / e weakly converges in L 2 ... P {{0,1) x R x M + ) to a 
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Similar arguments can be used in order to prove that up to a subsequence, (//) converges to a 
solution / of the Milne problem (|3.ip - (|3.2p - (|3.3p when I tends to +00. Indeed, if f\ admits the 
decomposition 

fi(x,p) = (ai(x)(p 2 + n) + bi(x)p x )(l + P) + wi(x,p), 



with 



J p x WiPdp = J (p 2 + n)wiPdp = 0, 



then the sequence (wi) is bounded in L 2 ... P (R + x R x R + ) and pointwise in x as in (|3.23j) . 

Pr(l+|p|J j-+P 

And so, up to a subsequence, {w{) converges weakly in L 2 3 P (R + x R x R + ) and also weak 

Pi-(l+|p|) T+T 

star in x, weak in p in L°°(R + ;Lp ,s 3 (R x R + )). The sequences (a/) and (6/) satisfy 

(jai + J p 2 x wiPdp^j =0, (jk + J p x (p 2 + n)wiPdp^j =0, 
so that 

jai(x) = - J p 2 x wi(x,p)Pdp + 7a/(0) + y p 2 x wi(0,p)Pdp, 

7&/(x) = - / ^(p 2 + n)wi(x,p)Pdp + 76/(0) + / p x (p 2 + n)wi(0,p)Pdp. 



It follows that the sequences (a/) and (bi) are uniformly bounded on R + , and so, up to a subsequence, 
converge weak star in x. The limit of (fi) is a weak solution to the problem. This weak solution 
belongs to D. 

We can now prove that b^ = 0. For this we notice that the discussion of this section up to (3.11) 
included, also holds for //, W being nonnegative on [0, /] because it is non increasing and vanishes at 
I. The discussion from (3.12) leading up to (3.15) is valid as well. But for // it holds that 6/(Z) = 0, 
and so (I3.14p taken at x = I leads to |6/oo| < ce~ vl . 
Take /J > a > 0. Using (|3,14p again implies that for all I > f3, 

Mx)\ < \b t (x) - b ioo \ + ce"^ < 2ce~ lia , x > a. 
It follows that 

|6(ar)| < 2ce- T > a , x > a. 
Hence 

lim b(x) = = b^. 



The uniqueness of the solution of the Milne problem (I3.ip - (l3.2p - (13.3p can be proven as follows. Let 
/ G D be solution to the Milne problem ()3.1[) - f)3.2|) - (j3.3j) with zero indatum at x = and zero 
energy flow. Let 

f(x,p) = a{x){p 2 + n)(l + P) + b{x)p x {l + P) + w(x,p) 

be its orthogonal decomposition. By (I3.1ip 

lim [ Px f\x,p)— ^- 5 dp = 0. (3.25) 
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Multiply the equation 

PxdJ = Lf, (3.26) 

by fj^pi integrate over ]0, +oo[xR 3 and use the spectral inequality. Then, 

If ~ P f +OQ f P 1 f ~ P 

77/ |Px|/ 2 (0,p) i : v dp+v / / w 2 (x,p) dpdx < -- lim / Pxf 2 (x,p) dp 

2 y Pa: <o 1 + p Jo J 1 + p 2 *-»+°° y 1 + p 

= 0. 

And so, 

/(0, -)=0, «;(-,•) =0. 
Equation f|3.26j) reduces to 

^7 = 0, 

so that together with /(0, •) = 0, it holds that a(-) = = 0. Hence / is identically zero. 
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